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Abstract 

We discuss the BLM scale fixing procedure in exclusive electroproduction processes in the 
Bjorken regime with rather large xb. We show that in the case of vector meson production 
dominated in this case by quark exchange the usual way to apply the BLM method fails due to 
singularities present in equations fixing the BLM scale. We argue that the BLM scale should 
be extracted from the squared amplitudes which are directly related to observables. 



Introduction 



The investigation of the quark-gluon dynamics through perturbative calculations is most useful to extract 
from experimentally measurable observables quantities such as parton distributions, generalized parton 
distributions (GPDs) and (generalized) distribution amplitudes (GDAs, DAs). Factorization theorems 
allow to calculate scattering amplitudes in a perturbative way, provided a renormalization scale and a 
factorization scale are choosen. Whereas the observables in the extensively studied inclusive reactions 
are in general related to an amplitude (as the case of inclusive DIS expressed as the imaginary part 
of the forward virtual Compton scattering reaction), exclusive cross sections which have been much 
studied recently are based on a factorization theorem at the amplitude level, and thus require to square 
the amplitude given by its perturbative expansion. Renormalization scale fixing has been the subject 
of intense studies and different strategies ^ |2] have been put forward to maximize the predictivity of 
theoretical studies through ensuring the smallness of corrections related to higher orders terms in the 
perturbative series. The phenomenological success of these proposals is quite impressive in a number of 
cases, mostly related to inclusive cross sections or jet physics. With the advent of next to leading order 
results, it has been advocated to use these procedures also in hard exclusive processes 011]. Even at 
the Born level, the hard meson electroproduction amplitude contains as. Therefore the choice of the 
renormalization scale ^xr is very crucial for the practical estimations of the observables related to meson 
electroproduction. The first study of a hard electroproduction amplitude including the analysis of the 
next-to-lcading orders (NLO) has been implemented by Belitsky and MuUer P] for tt"*" production, i.e. 
for 7*p TT+n process. This study used an appropriate continuation of the NLO calculations known for 
the electromagnetic pion form factor onto the case of the meson electroproduction process. 
In this work, focusing on the Brodsky-Lepage-Mackenzie (BLM) procedure P] ^, we examine in detail 
the consequences of the fact that exclusive processes considered in the regime where the quark GPDs 
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are dominant are factorized at the amplitude level and that the meson electroproduction amplitude is 
a complex function. Consequently, we are forced to apply the BLM procedure to the real and to the 
imaginary part of the scattering amplitude separately, which in general leads to two different scales. 
Moreover, we show that such a way of scale fixing, as has been done in for the tt— meson production, 
leads to unphysical results in case of vector meson production. We propose a way to modify the BLM 
procedure in order to avoid such difficulties. 



Basics of the BLM procedure 



The QCD factorization theorem [S] states that the amplitude of hard meson electroproductions can be 
written as 



i 1 

A = J dz J dx<^Aiiz,fJ,p)H{x,z,Q'^,fi%,fiji)F{x,^i%) = <^m'»H<»F, 



(1) 



where the parameters and are the factorization and renormalization scales, respectively. The 
scales /i^ and are in principle independent but often it is argued that they can coincide, fj,j^ — /ip. 
The arguments in favour of such assumption are discussed in, e.g. 0, and we adopt this also in the 
present paper (to simplify notation we omit below subscripts R and F). In Eq.l^, H is the hard part of 
amplitude which is controlled by perturbative QCD. The meson distribution amplitude describes the 
transition from the partons to the meson, and F denotes the GPDs which are related to nonperturbative 
matrix elements of bilocal operators between different hadronic states. 

The product (E) H (g) F in ^ is, generally speaking, independent of the particular choice of the 
parameter fi^. However, this independence is broken once we limit ourselves to the first few terms in 
an expansion over the coupling constant as- In this case, the theoretical ambiguity of the choice of 
the parameter /i^ emerges. The goal is to choose the parameter fi^ such as to ensure that pretty small 
contributions will arise from the next order corrections. Out of the several possible ways to hope to reach 
that goal, the BLM procedure 2 begins with separating out the terms which are proportional to the 
one-loop /3— function, Pq — 11 — 2/3Np, appearing in the NLO terms. The amplitude including the 
NLO corrections with separated terms proportional to (3o reads 
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y^LO(Q2)+iNLO,(^)(Q2 



+ 



(2) 



where the ellipsis stand for the terms of the NLO corrections which do not explicitly contain (Hq. In (|2J|, 
the value of the constant C depends on the kind of produced mesons. As pointed out in |S] the exact 
expressions, in the quark sector, with the NLO corrections may be obtained by a suitable substitution 
from the well-known results for the pion electromagnetic form factor. 
Due to the renormalization group equations the coupling constant takes the form 
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(/3o/4^)ln(AiVA^CD) 



(3) 



We insert this expression into the amplitude Q and then expand it in powers of asif^D- Retaining the 
terms which are proportional to /3o, we get 
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^LO(Q2)+iNLO,(/3)(Q2 
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The BLM procedure consists in the choice of such /i§ for which the whole term proportional to f3o in 
Eq. Q vanishes, i.e. 



2 / 2\ 00 



C-ln-^ 



(5) 



from which it follows that the BLM scale IJ^^lm equal to 
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_4NLO,(/3) 



(6) 



If the scattering amplitude is real, as in the case of the spacelike pion form factor, the above procedure 
leads just to the one ^j^blm scale. But already in the 7r+— meson electroproduction the scattering ampli- 
tude is a complex function and if one applies the BLM procedure separately both for real and imaginary 
parts "'j this results in two different scales. The situation starts to be even worse in the case of the vector 
meson electroproduction which we discuss in the next section. 



Extraction of the BLM scale from the ampUtudes 



We now focus on the numerical estimation of the renormalization scales extracted them from the ampli- 
tudes of the vector mesons electroproduction. We consider only spin non-flip quark GPD H of the nucleon 
target and neglect spin-flip GPD E. The reason is that: a) E function gives very small contribution to 
the scattering amplitude, b) its form is very model dependent. 
Consider the NLO terms of amplitude |5J containing the /3o coefficient 
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In Eqn. Q, the distribution amplitudes '^piz) and correspond to the p and hybrid mesons, 

respectively. The hybrid meson is a charge conjugation even state (J^*^ = 1 as studied in 8 . The 
function HP(x,£,,tmin) stand for the corresponding GPD's and is defined as 



HP{x, = ^ ( 1 1) - edH^x, ^, t) 



(8) 



The meson distribution amplitudes <I>p and may be understood as the asymptotic functions which 
are (see, for instance i8|): 



$p(z) = 6zz, 30zz(l - 2z). 



(9) 



Note that the NLO evolution effects for the meson distribution amplitudes seem to be small and we omit 
the consideration of such effects. 

The GPD's, in iQ, can be modelled using the Radyushkin model ;9 which ensures the agreement with 
the forward limit and the corresponding sum rules for the moments. According to this ansatz, the GPD's 

*The exchange of the on-shell (Hght-hke) gluon is entirely responsible for the imaginary part of the amplitude. This, 
however, does not break the factorization. 
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are expressed with the help of double distributions F'^{x, y; t): 



= j dyF\x+,y,t)~^-^^^ J dy F^{x^,y,t), (10) 



where 

For the double distribution F'^(X, Y;t), we assume the ansatz suggested by Radyushkin 

where the forward (anti)quark distribution is taken from the parameretization of ^Hl- Note that, in (|12|l . 
tmin is different from zero and is equal to — 4m^^^/(l — (see, for instance, JT]). A similar expression 
gives the anti-quark contribution. 

As shown in [12]^ , the definition of the double distribution is not completely compatible with the structure 
of the corresponding matrix elements; introducing D-tcrms restores the self-consistency of this represen- 
tation. Taking into account these D-terms with a factorized ^-dependence as in Eq. p2|l . the GPD's H1Q|I 
are modified into : 

HUx,U) - H''{x,tt)+e{^~\x\)^^^^, (13) 
where D{x/^,0) is given by (Ci?^^ are Gegenbauer polynomials) 

D{a) = -4(1 - a^)I^Cf^^{a) + 0.3Cl^^{a) + 0.icf\a)^ (14) 

with 

a=|, D{a)^-D(-a). (15) 

These D-terms do contribute to the interval — ^ < x < ^ of the GPD's. Besides, due to the anti-symmetric 
properties of (|14|l . the D-terms are important only for the charge conjugation odd vector meson {e.g. p) 
production amplitude rather than for exotic hybrid meson [J^'-' —1 ^) production amplitude [H]. 
As aforementioned, the amplitudes of the mesons electroproductions contain real and imaginary parts. So, 
we now come to the estimation of BLM scales for the real and imaginary parts of amplitudes, separately. 
Repeating the procedure pointed out in the preceding section, we derive that the BLM scales read 



where 



A^^.) = Q^ (16) 

^?/™)=e-/-(«)Q^ (17) 



fP (£) - 1^ - (18) 
f (.) - _ ^!tM. (19) 
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Figure 1: The p meson production: the BLM scale for the real part of the amplitude. 



for the p meson production, and 



(0' 



(20) 
(21) 



for the hybrid meson production. The explicit expressions for functions i?/^'^(^), -ff/™(^), Gf '^(^) and 
Gf™(^) in lO-lEU can be found in the Appendix. 

The investigation of the p meson scale with (|18ll shows that the extraction of the BLM scale from the 
expression for the amplitude meet difhculties. Indeed, as one can see on Fig. ^ the p meson function 
/fle(0 unphysical singularity owing to the fact that the denominator in p8l) may vanish (see 

the exact expression for iJ/^'^(f) in Appendix). Indeed, let us dwell on the equation (|28ll from the 
Appendix which is the denominator of Eq. p8|) . The integrand of (|28|l is a sign-changing function: the 
integrand is negative in the region — ^ < x < while it is positive in the regions — 1 < x < — ^ and 
^ < a; < 1. Qualitatively, it is clear that at some value of ^ the positive contribution to the whole integral 
(which is taken in the sense of the Cauchy pincipal value) will be equilibrated by the negative one. The 
numerical calculations show that the dangerous singularity appears for ^ « 0.27. This value depends 
on the parameterization of GPDs but the existence of a singularity is a model-independent result of our 
analysis. Concerning the function fj^{C) in (|19|l . as one can see from Fig. [21 this function is always 
analytical. 

For the hybrid meson scale, the situation is analogous. In this case, the sign of the integrand of Eq. H32|l 
of the Appendix is opposite to the one of the integrand in (|28|) : the integrand is positive in the region 
— f < X < ^ and is negative in the regions ^ < x < 1. Thus, we may expect that the whole integral H32|l 
can add up to the zeroth value. The BLM scale corresponding to the imaginary part of the hybrid meson 
production amplitude is again an analytical function. 

It is instructive to compare the BLM scale fixing for vector meson production with the one for the 
7r+ meson production 3 . In this second case such a singularity does not appear and equations fixing 
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Figure 2: The p meson production: the BLM scale for the imaginary part of the amphtude. 

BLM scale, both for real and imaginary parts of the scattering amplitude, are analytical. Indeed, the 
corresponding integral determining the BLM scale for the real part of the tt"*" meson production amplitude, 
i.e. 



ed 



(22) 



will never be equal to zero. 

Summarizing this section, the physical causes for the appearance of the singularity are the C-parity 
conservation and the factorization in hard reactions and the mathematical evidence is the sign-changing 
integrands of H28|l and H32() . Moreover, the vanishing of the first order term in the real part only does 
not imply that the scale of the gluon propagator vanishes. We thus see that the usual BLM scale fixing 
needs to be modified in the case of vector meson production. We propose to extract this scale directly 
from an observable, i.e. starting from the square of the scattering amplitude. 



Extraction of the BLM scale from the cross section 



In this section, we will now extract the corresponding BLM scales working with the squared amplitudes 
or, in other words, with the cross section. In this case, the BLM equation will be rewritten in the following 
form: 

4(^2) ^ |_^LO|2 _ 4(^2) ^jsRe^LO gffe^NLO,(/3) ^ cj^^LO cj^^NLO,(/3) | ^ 

(23) 
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Figure 3: The p meson production (quark contribution only): the BLM scales extracted from the squared 
amplitude. 



6,4- 



f^(^)=9/2 - [ G/^(^) G;'^^) + G;"^© G;"© ] / 
[ (G^^W + iz' (G;^^))' ] 



6,0- 



5,8- 



5,6- 
5,4- 
5,2- 
5,0- 



— I 

1,0 



0,0 0,2 0,4 0,6 0,8 

skewedness parameter ^ 



Figure 4: The hybrid meson production: the BLM scales extracted from the squared amplitude. 



7 



where we introduced the notations 

^NLO,(/3) 



C-ln 



21 



(24) 



From (|23|l . we can obtain for the p meson function /''(^) and for the hybrid meson function f^{£,) the 
following expressions: 



and 



^(0 



(25) 



(26) 



In (ESI and (EH, the structure functions "(C), ^,^^=(0 and Gf'"(0 are the same as they were 

defined for the BLM scales (jTH)) - (EJl- The curves of (ESJ and (HHJ are shown on Fig. |31and|3| 



Conclusions 



We have shown that the usual way of applying the BLM method to the scattering amplitude of exclusive 
vector meson electroproduction leads to equations fixing the BLM scale which are singular; we believe 
that this invalidates their straightforward use. The reason that this problem has not been taken care of 
in previous studies ^31 comes from the fact that these studies mostly concentrated on the simpler case of 
the meson form factor and of amplitudes which were rewritten in terms of this quantity. Let us stress that 
the singular behaviour of the equations is not specific to QCD. Indeed a simple calculation shows that the 
same behaviour would be obtained in a QED calculation where gluons are replaced by photons, quarks by 
electrons, and the p meson by a positronium bound state described by a conveniently defined distribution 
amplitude. We know that it is often advocated that the BLM procedure is easier to understand in the 
abelian case than in the non-abelian one, but its application to an exclusive process where the amplitude 
contains both a real and imaginary parts nevertheless suffers from the problem outlined in this paper. We 
have demonstrated that such singularities do not appear if the BLM procedure is applied to the square 
of the scattering amplitude, which is a quantity more closely related to an observable, the cross section, 
rather than the scattering amplitude itself. The phenomenological consequences for the vector meson 
electroproduction based on this new way of fixing of the BLM scale is studied in JI] . 
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Appendix: Typical functions for the determination of BLM scales 



The typical functions in terms of which the renormahzation scales corresponding to the p and hybrid 
mesons are rewritten read 
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for the hybrid meson production. 
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